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基于 引导 滤 疲 的 多 图 谱 医 学 图 像 分 割 
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摘要 :目的 为 了 有 效 的 利用 图 谱 的 先 验 信息 和 待 分 割 图 像 的 灰 度 信息 ,并 在 融合 标号 图 像 的 过 程 中 校正 配 准 引起 的 误差 ,得 到 


光滑 ,准确 的 分 割 结果 ,提出 了 一 种 新 的 基于 引导 滤波 的 多 图 


谱 医 学 图 像 分 割 方法 。 方 法 本 文 将 多 图 谱 配 准 与 引导 滤波 相 结 


合 。 该 方法 包含 4 个 部 分 :第 一 部 分 为 多 图 谱 配 准 ,通过 配 准 将 图 谱 中 存储 的 形状 先 验 信息 映射 到 待 分 割 图 像 ; 第 二 部 为 标号 融 
合 , 利 用 配 准 的 相似 性 作为 权重 ,将 形变 后 的 标号 图 像 融合 在 一 起 ;第 三 部 分 为 引导 滤波 ,利用 引导 滤波 引入 竺 分 割 图 像 的 灰 度 
信息 ,可 以 校正 配 准 引起 的 误差 ;最 后 通过 阔 值 处 理 , 得 到 最 终 的 分 割 结果 。 结 果 对 15 例 脑 部 MR 图 像 数据 中 的 海马 体 进行 分 


制 实验 ,左右 海马 体 分 别 达到 了 86% 及 87.4% 的 分 割 精 度 , 与 传统 的 标号 融合 算法 相 比 ,平均 分 割 精度 提升 了 2.4%。 结 论 本 


文 方法 结合 多 配 谱 配 准 与 引导 滤波 的 优势 ,提高 了 海马 的 分 割 精度 ,并 得 到 光滑 有 效 的 分 割 精 度 。 
关键 词 : 图 像 分 割 ;引导 滤波 ;多 图 谱 配 准 ;图 谱 先 验 ;标号 融合 ;海马 体 


Medical image segmentation based on guided filtering and multi-atlas 
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Abstract: A novel medical automatic image segmentation strategy based on guided filtering and multi-atlas is proposed to 
achieve accurate, smooth, robust, and reliable segmentation. This framework consists of 4 elements: the multi-atlas 


registration, which uses the atlas prior information; the label fusion, in which the similarity measure of the registration is used 
as the weight to fuse the warped label; the guided filtering, which uses the local information of the target image to correct the 


registration errors; and the threshold approaches used to obtain the segment result. The experimental results showed part 
among the 15 brain MRI images used to segment the hippocampus region, the proposed method achieved a median Dice 
coefficient of 86% on the left hippocampus and 87.4% on the right hippocampus. Compared with the traditional label fusion 
algorithm, the proposed algorithm outperforms the common brain image segmentation methods with a good efficiency and 


accuracy. 
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丘脑 海马 体 . 伏 隅 核 . 尾 状 核 等 关键 脑 结构 的 位 
置 ,体积 ,形态 等 的 变化 与 多 种 脑 部 疾病 息息相关 ,对 其 
精准 分 割 是 进行 定量 分 析 的 前 提 "。 脑 部 磁 共 振 图 像 
中 不 同 的 脑 组 织 相互 混 迭 ,结构 复杂 ,同时 ,由 于 磁 共振 
图 像 受到 伪 影 磁场 的 不 均匀 性 及 噪声 的 影响 ,使 得 在 
磁 共 振 岁 像 中 这 些 结构 边缘 模糊 .边界 难以 界定 ,因此 
PERRI AERA KIEKIS ,水平 集 的 方法 
等 多 种 方法 由 于 各 种 原因 难以 准确 分 制 ,例如 浆 值 法 
只 采用 了 图 像 中 像素 的 灰 度 信息 ,没有 使 用 像素 的 空间 
言 息 , 因 此 对 噪声 和 图 像 的 非 均 匀 性 敏感 。 区 域 生长 法 
对 每 一 个 需要 提取 的 目标 都 必须 手动 指定 一 个 种 子 点 ， 
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并 且 由 于 受 噪声 和 偏 移 场 的 影响 ,会 导致 区 域内 有 空洞 
其 至 是 根本 不 连续 的 分 制 结果 。 而 水 平 集 方法 “需要 
通过 周期 性 地 初始 化 水 平 集 孔 数 ,来 获得 稳定 的 曲线 演 
化 ,而 初始 化 水 平 集 函 数 计算 量 大 ,并 旦 很 难 确定 何 时 
和 如 何 重 新 初始 化 水 平 集 函 数 。 

而 基于 图 谱 配 准 的 分 割 方法 利用 手动 分 割 精 度 高 
的 优势 ,将 分 市 问 题 转化 为 配 准 问题 。Baillard 等 利用 
基于 单 图 谱 的 分 制 方法 分 割 出 脑 部 结构 。Heckemann 
等 "利用 29 组 图 谱 分 割 出 了 67 个 脑 部 结构 。Klein 等 了 
应 用 弹性 配 准 的 方式 分 割 3D MR 岁 像 中 的 前 列 腺 。 
Isgum 等 "使 用 前 向 搜索 法 来 减少 图 谱 数 以 提高 分 割 速 
度 。Jia 等 提出 把 基于 迭代 多 图 谱 方法 和 基于 树 形 结 
构 的 方法 相 结 合 以 充分 考虑 图 谱 间 及 待 分 割 图 像 间 相 
互 关 系 , 从 而 提高 了 图 像 的 分 割 精 确 。 以 上 基于 图 像 配 
准 的 分 割 方法 各 有 优势 ,但 得 到 的 结果 表面 不 光滑 ,在 
边界 处 会 出 现 阶梯 状 , 离 临床 应 用 还 存在 一 定 差距 。 本 
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文 将 多 图 谱 配 准 与 引导 滤波 相 结合 ,提出 了 一 种 基于 引 
导 滤 波 的 多 图 谱 医学 图 像 分 割 方法 。 


1 材料 与 方法 
1.1 实验 数据 

为 了 便于 和 有 目前 国际 上 主流 的 分 割 方法 进行 比较 ， 
我 们 采用 MICCAI 2013 竞 赛 中 的 公开 脑 部 数据 进行 了 
左右 海马 的 分 割 实验 , 该 数据 的 训练 集中 共有 35 组 脑 
部 Tl 加 权 磁 共振 三 维 数据 ,包含 灰 度 图 像 及 对 应 的 专 
家 手动 勾画 的 标号 图 像 。 我 们 随机 的 选择 20 组 作为 图 
谱 , 其 余 15 组 作为 测试 。 计 算 重 琶 率 Dice Fil Hausdorff 
距离 用 于 验证 分 制 方法 的 精度 ,其 定义 分 别 为 : 

2V(A NB) 

V(A)+ V(B) 

HD=max(Hi, H2) 

其 中 分 别 表示 自动 分 割 结 果 和 专家 手动 勾画 结果 
的 体积 ,V(ANB ) 表示 自动 分 割 结果 和 专家 手动 勾画 结 
果 重 登 部 分 的 体积 ;Hi=maxsa(minsesd(a,b)) ,Hs=maxvs 
(minsead(a,b)),d(a,b) 表 示 a,b 两 个 点 的 欧 氏 距离 。 
1.2 算法 流程 

本 文 方法 由 图 像 预 处 理 .多 图 谱 配 准 标号 融合 . 引 
导 滤 波 姜 值 处 理 五 部 分 组 成 。 本 文 方法 系统 流程 网 如 
图 1 所 示 。 


Dice(A, B)= 


Training Images | 


ilas Probability] 


Guided Filter 
Threshoalding 
Final Result 
图 1 系统 流程 图 


Fig.1 Flowchart of the proposed method. 


1.2.1 图 像 预 处 理 图 像 预 处 理 是 图 像 分 割 的 第 一 步 ， 
在 前 期 的 工作 基础 中 ,我 们 发 现 脑 元 和 偏 移 场 会 对 图 像 
的 分 割 有 一 定 的 影响 ,我 们 用 BET 算 法 "去 除 脑 却 ; 用 
N4 算 法 "去 除 磁 共振 图 像 中 的 偏 移 场 ;进行 灰 度 归 一 
化 完成 图 像 的 预 处 理 。 下 图 2 为 预 处 理 后 的 脑 部 磁 共 
振 图 像 。 


图 2 预 处 理 后 的 脑 部 磁 共振 图 像 
Fig.2 Preprocessed brain MR image. 


1.2.2 配 准 在 基于 图 谱 配 准 的 分 割 方法 中 ,图 谱 的 选 
择 对 分 割 结果 有 一 定 的 影响 ,为 了 降低 图 谱 选 择 的 不 确 
定性 ,我 们 选取 多 张 图 谱 , 图 谱 由 两 部 分 组 成 :图 谱 灰 度 
图 像 1 和 影像 专家 手动 分 制 好 的 标号 图 像 L,L 是 一 个 
一 值 图 像 。 首 先 , 将 图 谱 灰 度 图 像 I 配 准 到 待 分 割 图 像 
T 得 到 变形 场 。 然 后 用 变形 场 对 图 谱 标号 图 像 L 变 形 ， 
得 到 变形 后 的 标号 图 像 L', 最 后 融合 变形 后 的 标号 图 像 
LL; 即 可 得 到 最 终 的 分 割 | 结果 。 

1.2.3 概率 图 谱 传统 方法 在 融合 时 一 般 需 要 重新 计算 
变形 后 的 图 谱 与 待 分 制图 像 的 相似 性 ,如 互信 息 量 、 互 
相关 系数 . 归 一 化 的 互相 关 或 者 其 它 的 相似 性 测度 ,这 
种 方法 额外 增加 了 计算 量 , 比较 耗 时 ,而 且 这 种 整体 的 
相似 性 测度 并 不 能 很 好 的 度量 每 个 像素 点 在 配 准 、 分 割 
时 的 作用 。ANTS 配 准 方法 利用 局 部 相似 性 作为 配 准 
的 相似 性 测度 ,并 且 对 变形 场 施加 了 有 效 的 约束 及 平滑 
可 以 达到 比较 好 的 配 准 效 果 "。 我 们 以 ANTS 配 准 的 
局 部 相似 性 作为 融合 时 的 权重 ,相似 性 测度 高 的 像素 
点 ,说 明 配 准 更 加 准确 ,对 最 后 的 分 割 结果 影 响 就 越 
大 。 局 部 权重 概率 图 谱 的 计算 公式 为 : 


id i L'; 
Le Din XE (1) 
Se 
其 中 wo 为 ANTS 配 准 的 局 部 相似 性 测度 ,为 变形 
后 的 标号 图 像 。 


1.2.4 引导 滤波 引导 滤波 ”于 2013 年 提出 , 它 是 一 种 
边缘 保持 滤波 ,可 以 实现 图 像 边 缘 的 平滑 细节 描 强 A 
及 图 像 融合 去 噪 等 功能 。 主 要 原理 是 通过 一 幅 引导 图 
像 对 输入 图 像 进 行 滤波 处 理 ,输出 图 像 在 保留 输入 图 像 
整体 特征 的 时 候 还 能 充分 获得 引导 图 像 的 边缘 特征 。 
该 算法 涉及 一 个 引导 图 像 1, 输 入 图 像 P, 输 出 图 像 O。 

下 图 显示 了 引导 滤波 的 结果 ,图 3A 为 133x138 大 
小 的 脑 部 MR 图 像 。 将 图 3A 作为 引导 滤波 的 引导 图 像 
和 输入 图 像 ,图 3B 为 滤波 输出 图 像 。 对 比 图 3A 和 图 
3B 可 以 看 出 ,引导 滤波 会 保留 图 像 的 边缘 特征 ,但 会 弱 
化 图 像 的 细节 。 改 变 输入 图 像 , 图 A 仍 为 引导 图 像 ,将 
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图 3C 作 为 输入 图 像 ,而 图 3D 为 引导 滤波 后 的 结果 。 观 
察 图 3C, 可 以 看 出 图 3C 受 噪声 影响 ,并 且 在 图 像 边 缘 处 
有 孤立 点 ,对 比 图 3C 图 3D 可 以 看 出 ,引导 滤波 后 ,图 像 


| 


k 


m u 
S> 


图 3 引导 滤波 的 结果 


的 噪声 被 消除 ,并 且 那 些 孤立 的 点 和 图 像 边缘 连接 到 一 
起 。 如 前 所 述 ,输出 图 像 在 保留 输入 图 像 整 体 特征 的 时 
候 还 能 充分 绑 得 引导 图 像 的 边缘 特征 。 


Fig.3 Two examples of guided filtering. A, C: Two input images; B, D: The output images of the former two images independently. 


1.2.5 BARR 在 本 文中 ,由 公式 (1) 得 到 局 部 权重 概 
率 图 谱 , 将 得 到 的 概率 图 谱 工 进行 引导 滤波 ,引导 图 像 
为 待 分 割 图 像 I。 通 过 这 种 局 部 线性 模型 ,我们 将 引导 
图 像 ( 待 分 割 图 像 ) 中 的 结构 信息 引入 到 输出 图 像 ,滤波 后 
的 概率 图 谱 ( 输 出 图 像 ) 记 作 FL, 则 最 终 的 分 割 结果 了 为 : 


1,FL>M 
R=\0,FL<M (5) 


2 结果 与 讨论 

为 验证 本 文 方法 的 有 效 性 ,我 们 与 目前 国际 上 主流 
的 几 种 方法 进行 了 比较 ,分 别 为 STAPLE™, Spatial 
STAPLE "7 , Major Voting "9 , Weight Voting "和 
SIMPLE", E S X Windows 7 Service Pack 1, 
Intel(R) Core(TM) i5-3330 CPU,8.0 G 内 存 。 

Fe Arie OR : R 1 A S ARAE 
对 左 Ar ie Fh oP BSH AR AY FAK. STAPLE, Spatial 
STAPLE, Major Voting, Weight Voting 和 SIMPLE 五 种 
方法 对 左右 海马 分 制 的 重 闭 率 分 别 为 0.84 BQ 0.85 FE 
右 ,而 本 文 方法 可 以 达到 0.86 及 0.87, 提 高 了 2.4%。 图 
4 给 出 了 重 登 率 盒 状 图 (第 一 行 ) 和 Hausdorff 距 离 盒 状 
图 (第 二 行 ), 红 线 为 中 位 数 。 从 盒 状 图 中 可 以 看 出 本 文 
方法 重 登 率 和 距离 值 都 比较 集中 ,而 其 它 方法 的 结果 则 
比较 分 散 , 说 明 本 文 方法 具有 较 好 的 一 至 性。 重 装 率 越 
高 说 明 分 制 精 度 越 高 ,本 文 方法 的 重 闭 率 中 位 数 在 盒 状 
图 的 偏 上 方 ,说 明 对 于 大 部 分 测试 图 像 都 可 以 达到 较 好 
的 分 割 。Hausdorff 距 离 值 越 小 说 明 分 割 精度 越 高 ,本 


表 1 不 同方 法 左 \ 右 海马 分 割 结果 与 手动 分 割 结果 的 重 垒 率 
Tab.1 Dice similarity for the various segmetation as compared 
with manual segmentations 


Method Left Hippocampus Right Hippocampus 
STAPLE 0.83+0.07 0.84+0.06 
Spatial STAPLE 0.84+0.06 0.84+0.06 
Major Voting 0.83+0.07 0.84+0.06 
Weight Voting 0.84+0.06 0.85+0.05 
SIMPLE 0.84+0.06 0.85+0.05 
Our 0.86+0.04 0.87+0.03 


图 SA~F 分 别 为 本 文 方法 及 STAPLE, Spatial 
STAPLE, Major Voting, Weight Voting, SIMPLE 对 左 、 
右 海马 分 割 的 对 编号 为 1006 的 数据 的 分 制 结 果 。 绿 色 
为 专家 手动 分 割 结果 ,红色 为 不 同方 法 的 分 制 结果 。 从 
A 中 可 以 看 出 本 文 方法 与 专家 手动 分 割 结果 更 加 接近 ， 
特别 在 脑室 附近 (黄色 箭头 所 示 ) ,其实 方法 都 出 现 了 错 
误 的 分 制 ,而 本 文 方法 在 通过 引导 滤波 引入 竺 分 割 图 像 
的 信息 ,成 功 的 将 演化 曲线 吸引 到 了 真实 的 海马 边界 。 

比较 各 方法 三 维 重 建 的 结果 (如 图 6) ,用 其 它 方法 
所 得 到 的 结果 存在 明显 的 阶梯 状 效果 ,在 边界 出 会 出 现 
一 些 毛刺 或 凹陷 现象 ,而 新 方法 利用 引导 滤波 可 以 得 到 
光滑 .连续 的 边界 ,与 脑 部 结构 的 真实 情况 更 加 接近 , 便 
于 医生 观察 体积 与 形态 的 变化 。 


3 结论 


大 脑 内 部 结构 的 位 置 体积 与 形态 变化 与 多 种 疾病 


SOME Hausdorff HE S PACE REL i FA, E 
样 可 以 说 明 对 于 大 部 分 测试 图 像 本 文 方法 可 以 达到 较 
好 的 分 割 , 红 点 为 出 格 点 , 即 偏离 平均 值 较 远 的 点 ,本 文 
方法 没有 出 格 点 ,而 其 它 方法 都 存在 出 格 点 ,同样 说 明 
了 本 文 方法 具有 较 好 的 一 致 性 。 


息息相关 ,而 要 确定 和 分 析 这 些 变 化 首先 需要 对 这 些 脑 
结构 进行 精确 的 分 割 。 我 们 对 目前 脑 内 结构 分 割 方法 
进行 了 详细 分 析 ,并 针对 这 些 方法 存在 的 缺点 和 不 足 ， 
提出 了 一 种 新 的 基于 引导 滤波 的 多 图 谱 医 学 图 像 分 割 
方法 。 新 方法 中 包含 两 个 重要 部 分 :多 图 谱 配 准 和 引导 
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图 4 不 同方 法 分 割 结果 重 琶 率 (第 一 行 ) 和 Hausdorff 距 离 ( 第 二 行 ) 的 箱 线 图 

Fig.4 Box-plot of performance metric for segmentation methods, as compared against manual segmeatation . 
The top row contains the dice similarity, and the bottom row contains the Hausdorff distance. The left 
column is the left hippocampus result and the right column is the right hippocampus. 


5 不 同方 法 对 编号 为 1006 的 数据 的 分 割 结果 

Fig.5 Sagital view of the various segmetation results of MICCAI image. A-F: STAPLE, Spatial STAPLE, Major Voting, 
Weight Voting and SIMPLE. Automatic segmentation results are illustrated in red contour, the manually-delineated the 
hippocampus boundary shown by the green contour. 


滤波 。 概 率 图 谱 的 准确 性 将 直接 影响 到 分 制 的 精度 , 利 ”参考 文献 : 
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Fig.6 Surface rendering of ground truth and segmentation results. Three-dimensional hippocampus segmentations 
obtained by six methods. A is the manual segmentation; B-G: the proposed method, STAPLE, Spatial STAPLE, Major 


Voting, Weight Voting, SIMPLE. 


[3] Lu XQ, Wu JS, Ren XY, et al. The study and application of the 
improved region growing algorithm for liver segmentation [J]. 
Optik, 2014, 125(9): 2142-7. 

[4] Bai PR, Liu QY, Li L, et al. A novel region-based level set method 
initialized with mean shift clustering for automated medical image 
segmentation[J |. Comput Biol Med, 2013, 43(11): 1827-32. 

[5] Baillard C, Hellier P, Barillot C. Segmentation of brain 3D Mr 


images using level sets and dense registration[ J]. Med Image Anal, 
2001, 5(3): 185-94. 

[6] Heckemann RA, Hajnal JV, Aljabar P, et al. Automatic anatomical 
brain MRI segmentation combining label propagation and decision 
fusion| JJ. Neuroimage, 2006, 33(1): 115-26. 

[7] Klein S, Van Der Heide UA, Lips IM, et al. Automatic segmentation 


of the prostate in 3D Mr images by Atlas matching using localized 
mutual information| J]. Med Phys, 2008, 35(4): 1407-17. 


[8] Isgum I, Staring M, Rutten A, et al. Multi-Atlas-Based segmentation 
with local decision Fusion-Application to cardiac and aortic 
segmentation in CT scans|J|. IEEE Trans Med Imaging, 2009, 28 
(7): 1000-10. 

[9] Jia HJ, Yap PT, Shen DG. Iterative multi-atlas-based multi-image 
segmentation with tree-based registration|J il; Neuroimage, 2012, 59 
(1, SI): 422-30. 

[10] Popescu V, Battaglini M, Hoogstrate WS, et al. Optimizing parameter 

choice for FSL-Brain Extraction Tool (BET) on 3D T1 images in 

multiple sclerosis[J]. Neuroimage, 2012, 61 (4): 1484-94. 


[11] Tustison NJ, Avants BB, Cook PA, et al. N4ITK: improved N3 bias 


correction|J |. IEEE Trans Med Imaging, 2010, 29(6): 1310-20. 
[12] Avants BB, Tustison NJ, Song G, et al. A reproducible evaluation of 
ANTs similarity metric performance in brain image registration[J]. 
Neuroimage, 2011, 54(3): 2033-44. 
[13] He KM, Sun J, Tang XO. Guided image filtering [J]. IEEE Trans 
Pattern Anal Mach Intell, 2013, 35(6): 1397-409. 
[14] Warfield SK, Zou KH, Wells WM. Simultaneous truth and 


performance level estimation (STAPLE): An algorithm for the 
validation of image segmentation |J]. IEEE Trans Med Imaging, 
2004, 23(7): 903-21. 

[15] Commowick O, Akhondi-Asl A, Warfield SK. Estimating a reference 
standard segmentation with spatially varying performance 
parameters: local MAP STAPLE [J]. IEEE Trans Med Imaging, 
2012, 31(8): 1593-606. 

[16] Aljabar P, Heckemann RA, Hammers A, et al. Multi-atlas based 
segmentation of brain images: Atlas selection and its effect on 
accuracy|J |. Neuroimage, 2009, 46(3): 726-38. 

[17] Sjöberg C, Ahnesjs A. Multi-atlas based segmentation using 
probabilistic label fusion with adaptive weighting of image 
similarity measures [J]. Comput Methods Programs Biomed, 2013, 
110(3): 308-19. 

[18] Langerak TR, Van Der Heide UA, Kotte A, et al. Label fusion in 
atlas-based segmentation using a selective and iterative method for 
performance level estimation(SIMPLE) [J]. IEEE Trans Med 
Imaging, 2010, 29(12): 2000-8. 


(编辑 :经 We) 


